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INTRODUCTION 

't 

This  report  presents  the  performance  of  a new  modulation  technique 
which  has  been  called  mixed- base  modulation  (MBM).  This  technique 
provides  significant  performance  advantages  over  other  common  modu- 
lation techniques  for  applications  where  both  power  and  bandwidth  are 
limited.  Its  performance  for  bandspreading  factors  of  two  to  five 
relative  to  single  - sideband  is  reasonably  close  to  the  ultimate  Shannon 
bounds.  Most  available  techniques  are  very  poor  in  this  region. 

The  viewpoint  which  allowed  synthesis  of  this  technique  will  be 
given  before  the  basic  approach  and  its  performance  are  discussed.  A 
section  comparing  the  theoretical  performance  of  MBM  and  that  of 
other  common  techniques  is  followed  by  a discussion  of  some  of  the 
practical  considerations. 
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THE  TRANSMISSION  PROBLEM 

The  problem,  to  be  considered  is  illustrated  in  Figure  1,  An  input 
signal,  m(t),  is  to  be  transformed  in  some  manner  to  a second  signal, 
h(t),  to  which  independent  white  gaussian  noise,  n(t),  is  added  to  produce 
the  received  signal,  r(t).  This  signal  is  then  transformed  into  an  esti- 
mate of  m(t),  m(t).  We  will  take  the  liberty  of  calling  the  device  that 
transforms  m(t)  into  h(t)  the  modulator  and  the  device  which  trans- 
forms r(t)  into  A(t)  the  demodulator.  The  fidelity  criteria  will  be  the 
normalized  mean- squared  error,  e,  between  m(t)  and  m(t)  where  a 
known  time  delay,  Td,  is  allowed.  Thus: 

_ E[(m(t)  - m(t+Td))8] 

G E[ma(t)]  ( 

The  objective  will  be  to  choose  the  modulation  and  demodulation  pro- 
cesses to  produce  a small  value  of  c when  the  spectral  characteristics  of 
the  source,  and  the  transmitted  signal,  m(t)  and  h(t),  and  the  noise 
spectral  height,  N / 2,  are  specified. 

To  simplify  the  problem  we  will  assume  the  input  data  to  be  a 

stationary  process  with  an  ideal  flat  bandlimitted  spectrum  of  bandwidth, 

B . We  then  can  represent  the  process  with  time  samples  of  the  input 

taken  nt  intervals  of  T : 

s 


T = 1/2B  . 


If  the  output  of  the  modulator  is  also  a sequence  of  samples  at 
intervals,  T^,  they  can  be  represented  by  time  functions  such  that  the 


* The  objective  of  centering  the  spectrum  of  h(t)  for  propagation  will  be 

discussed  separately  in  Section  VHI. 
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output  or  channel  bandwidth  will  be  B = 1/2T  . Thus,  the  bandspreading, 

c o 

fi  = B /B  , will  be  the  ratio  of  sampling  rates  between  the  output  and  in- 
c v 

put.  Further,  if  the  receiver  first  passes  the  noisy  signal  through  an 

ideal  sharp  cutoff  filter  of  bandwidth  B , and  the  result  is  sampled  at 

c 

the  appropriate  times,  the  problem  reduces  to  a transformation  of 
independent  random  variables  as  indicated  in  Figure  2 where  i is  the 
time  index  and  the  variance  of  the  noise  random  variables,  n , is: 


E[n.s]  = B N = RN  B . 

1 CO  o V 


It  is  worth  noting  that  our  setup  of  the  problem  differs  from  that 
taken  by  Ford  [1]  which  models  the  input  source  as  a finite  dimensional 
Markoff  process  and  allows  control  theory  techniques  discussed  by 
Snyder  [2]  to  synthesize  the  optimum  modulator-demodulator  pair.  Ford 
was  unsuccessful  in  evaluating  the  optimum  non-linear  pair  and  discusses 
instead  some  of  the  properties  of  the  optimum  time-variant  linear  combi- 
nations. It  should  be  pointed  out  also  that  the  real-time  (that  is  = 0) 
evaluation  of  error  implicit  in  the  model  provided  by  Snyder  is  usually 
unrealistic  in  communication  problems. 

The  simplest  case  to  consider  is  that  of  unity  bandspreading  and  a 
memoryless  modulator.  This  implies  that  the  modulator  is  simply  an 
operator  that  modifies  the  value  of  m.  to  produce  h.  and  the  receiver 
estimates  nrv  from  each  input,  r..  Since  the  received  signals  are 
independent  from  sample  to  sample  the  receiver  cannot  make  use  of 
other  received  samples.  If  we  define: 


and 


[n\] 


m.  = 0. 


Cril  = 


Cni  + 


u 

0 

L 


o 


then  one  may  wonder  if  a good  choice  of  non-linear  operators  might 
significantly  reduce  e (as  compared  with  e for  linear  operators)  for  a 
specified  value  of  receiver  input  signal  to  noise  ratio,  SNR.  = 

P /N  B = 2E  /N  . 

S O V so 

where: 

P = E [ h3  ] 
s x 

E = P T . 
s s s 

In  Appendix  I an  integral  equation  is  derived  for  the  optimum  modu- 
lator-demodulator pair  which  depends  upon  the  probability  density 
function  of  the  input.  An  explicit  solution  was  obtained  only  for  a 
gaussian  density  function. 

Lechleider  [3,4]  and  Haddad  [5,6]  have  considered  the  trans- 
mission of  the  memoryless  time  discrete  variables  and  derived  some 
results  for  the  "small  noise"  case  using  a truncated  Taylor  series 
expansion  for  the  non-linear  operator.  Both  have  extended  the  analysis 
to  consider  the  memoryful  case  where  each  output  sample  is  a function 
of  more  than  one  input  sample. 

Although  the  case  of  unity  bandspreading  is  interesting,  it  can  be 
shown  by  information  theoretic  arguments  that  where  the  input,  m(t), 
is  gaussian,  minimum  mean-square  error  is  achieved  when  the  operators 
are  linear  and  memoryless.  Goblick  [ 7 ] and  Brown  and  Palermo  [ 8 1 
show  that  the  maximum  achievable  output  signal-to-noise  ratio,  SNR  , 
for  a gavxssian  source  is: 

P B /B  p 

SNR  = 1/e  = (l  + r—')  C V = f 1 + SNR. /p  ] . (2) 

o \ In  d / l 

o c 

When  the  bandspread  ratio,  B /B  , equals  one,  the  maximum  value  of 

c v 
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SNR  is  achieved  when  both  modulator  and  demodulator  are  linear  and 

o 

memoryless  and  a voltage  gain  of  SNR./(1  + SNR.)  is  used  in  the 
demodulator. 

If  the  source  is  non-gaussian  the  calculation  of  the  rate  distortion 
bound  becomes  more  difficult.  However,  by  use  of  bounds  given  by 
Shannon  and  Weaver  [9  ] we  may  determine  bounds  for  SNRq  to  be: 

(1  + SNR./g)8  sSNR  sK(l  + SNR./p)3  (3) 

1 O 1 

where: 


2-r  ^ lo§5  £(m)dm*l 

i_  ' _i 


(2we)  ' )m2f(m)dm 


and  as  before  f(m)  is  the  probability  density  function  of  m.  If  the  source 
density  function  is  uniform  rather  than  gaussian,  K = 1.42.  Thus,  the 
best  possible  modulator-demodulator  for  uniform  density  and  bandspread- 
ing of  one  can  at  most  gain  1.  5 dB  over  the  linear  combination  previously 
discussed. 

From  equation  (3)  we  can  conclude  that  to  supply  a value  of  SNR^  of 
40  dB  an  input  signal  to  noise  ratio,  SNR  , of  between  38.  5 and  40  dB 
is  required  for  the  case  of  a uniformly  distributed  source.  An  output 
signal-to-noise  ratio  of  40  dB  is  achieved  with  an  SNR  of  40  dB  by  simply 
transmitting  the  modulating  signal,  m(t),  and  using  r(t)  as  the  estimate. 
Thus,  the  ultimate  results  are  achievable  (or  nearly  so)  by  use  of  a linear 
memoryless  modulator-demodulator  pair  where  all  of  the  error  is  due  to 
link  error. 

This  leads  one  to  wonder  if  relatively  simple  schemes  might  provide 
good  results  for  larger  bandspread  ratios.  Equation  (3)  indicates  that  if 


r 


f 
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a bandspread  ratio  of  two  is  employed  a scheme  exists  for  providing  a 

SNR  of  40  dB  with  23  dB  SNR.  for  a uniform  source.  This  performance 
o 1 

is  very  significant  since  a standard  FM  system  requires  a 22  dB  SNR. 

for  40  dB  SNR  with  a bandspread  ratio  of  20  and  a PSK  system  requires 

° * 

21  dB  SNR  for  40  dB  SNR  for  a bandspreading  of  14.  Thus,  a choice 
i o 

of  operators  must  exist  that  provides  performance  close  to  that  of  FM 
and  PSK  with  a bandspread  ratio  of  1/10  or  1/7  of  the  respective  system. 

The  problem  of  synthesizing  a modulator-demodulator  pair  for  a g 
of  two  utilizing  a memoryless  operator  is  illustrated  in  Figure  3.  Each 
input  pulse  is  passed  through  two  memoryless  transformations,  0j  (m) 
and  0?  (m)  to  produce  two  output  pulses  hx  and  h?  for  each  input  pulse. 

The  received  signals  rx  and  rp  are  the  sum  of  these  values  and  independent 
noise.  The  minimum  mean- square  error  estimate  for  m,  m,  is: 

m = E(m  j rx  , r2)  = OglTj  , r?)  (4) 

Expressions  for  and  e are  given  in  Appendix  1 for  this  case.  Although 
in  principle  these  equations  can  be  used  to  provide  a solution  for  the 
optimum  choice  of  01(m),  0g(m)  and  0g(r1,  rg),  the  computational  diffi- 
culties appear  to  be  very  great  even  if  a computer  is  employed.  Thus 
the  selection  of  operators  for  higher  bandspreadings  than  one  requires  a 
different  approach. 


! 


* For  more  details  see  Appendix  5. 

**  The  subscript  i has  been  dropped  on  m,  h and  r in  the  memoryless 
case. 
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A notion  which  was  introduced  quite  early  in  the  information  theory 
literature  by  Shannon  [ 10]  (and  suggested  earlier  in  much  less  explicit 
fashion  by  Kotelnikov  [111)  and  has  been  discussed  more  recently  by 
Wozencraft  and  Jacobs  T1Z]  and  Sakriscn  [13]  is  that  of  viewing  the 
problem  geometrically. 

Let  us  consider  the  input,  m,  to  be  constrained  to  a finite  range 
(say  ± A)  and  assume  that  we  wish  to  choose  the  two  memoryless  operators, 
01(m)  and  03(m)  for  a bandspread  factor  of  two.  Any  choice  of  01(m)  and 
02(m)  can  be  viewed  as  a mapping  of  a line  representing  the  values  of  m 
onto  a plane  with  coordinates  hx  and  h3  as  illustrated  in  Figure  4.  Further 
the  pair  of  received  signals,  r and  rp  represent  a point  in  the  plane  that 
differs  from  the  transmitted  pair  of  signals  by  the  two  noise  components, 
nT  and  np  which  are  independent  and  gaussian.  The  task  of  the  demodulator, 
0 (t-j,  rp)  is  to  estimate  which  value  on  the  line  was  sent  given  a point  in 
the  plane.  Perhaps  the  simplest  rule  (but  not  necessarily  the  optimum  rule) 
for  the  receiver  is  to  select  the  point  on  the  line  closest  to  the  received 
point.  This  would  correspond  to  assuming  that  the  most  likely  noise 
vector  occurred  and  as  such  would  be  the  maximum  likelihood  estimate. 
Assuming  this  type  of  receiver,  how  should  we  draw  the  line? 

The  first  observation  we  can  make  is  that  if  the  noise  is  "small"  so 
that  we  do  not  confuse  one  major  segment  of  the  line  with  another  and  we 
choose  to  represent  the  noise  components  as  one  tangent  to  the  line  at  the 
estimated  value  and  one  perpendicular  at  the  estimated  value,  the  only 
error  will  be  that  due  to  the  component  tangent  to  the  line.  The  error 
in  the  estimate  of  m,  however,  depends  upon  how  different  the  estimated 
value,  m,  is  from  the  transmitted  signal  m^.  If  we  uniformly  mapped  m 
onto  a line  that  was  ZOA  in  length  (or  stretched  by  ten)  then  the  error  in 
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m normalized  to  the  full  scale  of  m will  be  one  tenth  the  size  of  the  noise 
component.  If,  however,  the  line  is  only  2A  in  length  the  error  in  m for 
"small"  noise  will  be  equal  to  the  noise.  Even  if  we  constrain  the  power 
in  the  two  pulses: 

Pg  = E[h3  + h3  ] = E[m3  ] (5) 

we  can  make  the  normalized  error  in  the  estimate  as  small  as  we  wish  for 

the  "small  noise"  case  by  making  the  length  of  the  line  sufficiently  long. 

In  this  case  the  gain  in  SNR^  for  a specified  SNR.  is  equal  to  the  square 

of  stretch  factor.  The  problem  with  this  is  that  for  a fixed  P we  cannot 

make  the  line  longer  without  reducing  the  distance  between  locations 

representing  largely  different  values  of  m.  When  the  noise  is  such  that 

errors  of  this  nature  (anomolies)  occur  frequently,  a thresholding  effect 

is  observed  in  the  SNR  versus  SNR  curve  and  the  closer  the  line  is 

o i 

packed  in  the  two  dimensional  space  the  smaller  the  noise  (or  the  larger 

SNR  ) necessary  to  cause  thresholding, 
i 

This  then  instructs  us  as  to  how  we  should  draw  the  line.  We  should 
make  the  line  sufficiently  long  to  provide  the  high- signal  to  noise  improve- 
ment that  we  wish,  but  make  it  such  that  the  probability  of  anomoly  is  as 
small  as  possible.  It  should  be  noted  that  the  arguments  presented  place 
no  particular  desirability  on  making  the  line  continuous.  This  is 
emphasized  since  references  [121  and  [13]  both  place  the  constraint  of 
continuity  upon  the  mapping.  Figure  5 shows  three  mappings  which  seem 
as  though  they  would  be  pretty  good.  Mappings  a and  b contain  no  dis- 
continuities whereas  mapping  c does  contain  discontinuities.  One  might 
expect  that  the  thresholding  effect  in  c would  be  sharper  than  in  a and  b 
where  an  anomoly  is  somewhat  fuzzy.  In  all  cases  a trade  between  high 
SNR  gain  and  threshold  can  be  made  by  putting  in  more  or  less  revolutions, 
sawteeth,  or  lines  as  the  case  may  be. 


10 
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If  a bandspread  ratio  greater  than  two  is  used  with  memoryless 
operators  the  mapping  problem  is  one  of  mapping  a line  into  more  than 
two  dimension  with  the  same  notions  of  stretching  and  spacing.  Timor 
[14  1 has  investigated  the  problem  of  mapping  a continuous  line  in  five 
dimensions  such  that  points  on  the  line  corresponding  to  values  of  m 
differing  by  more  than  a specified  value,  £m,  are  separated  by  a dis- 
tance ^r  or  greater.  Under  this  constraint  he  maximizes  the  stretch 
factor  and  estimates  resulting  performance. 

If  one  wishes  to  consider  the  use  of  memoryless  operators  then  the 
geometrical  object  to  be  mapped  is  no  longer  a line  but  a hyper- square. 
Thus  if  a bandspread  of  two  is  allowed  with  a memory  of  two  samples  the 
problem  is  one  of  mapping  a square  into  four  dimensions. 


1 
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MIXED  BASE  MODULATION 


By  use  of  the  geometric  approach  one  can  conjure  up  many  interest- 
ing configurations  for  mapping  lines  into  two  or  three  dimensions.  The 
type  of  mapping  shown  in  Figure  5c  is  particularly  interesting  since  it 
looks  as  if  it  should  be  good,  it  is  easily  implemented  and  its  perform- 
ance is  relatively  easy  to  evaluate.  Further,  it  can  be  easily  extended 
to  higher  dimensions  by  noting  that  the  location  of  parallel  lines  in  n 
dimensions  can  be  identified  by  the  location  of  points  in  n-1  dimensions. 
If  the  length  of  each  line  segment  is  2A  then  the  total  length  is  2A  times 
the  number  of  such  line  segments.  To  minimize  the  problem  of  anomaly 

(or  minimize  the  threshold  SNR  ) one  must  maximize  the  minimum  dis- 

1 

tance  between  these  points.  This,  however,  is  identical  to  the  usual 
discrete  signal  design  problem. 


The  implementation  of  such  a mapping  can  be  explained  by  consider- 
ing an  example  of  mapping  the  input  into  eight  lines  as  shown  in  Figure  6. 
The  two  operators,  (^(m)  and  0?(m)  required  to  achieve  this  mapping 
are  shown  in  Figure  7.  The  second  operator  can  be  recognized  as  an 
8-ary  quantizer  and  the  first  operator  as  one  that  provides  the  remainder 
or  quantizing  error  associated  with  the  8-ary  quantizer.  A simple  way  to 
synthesize  such  operators  is  to  use  a 3 bit  analog-to-digital  converter 
followed  by  a 3 bit  digital -to- analog  converter  to  produce  the  8-ary  pulse. 
The  remainder  may  be  obtained  by  subtracting  the  input  from  the  pulse 
output  as  shown  in  Figure  8 and  amplifying.  A technique  similar  to  this 
was  recently  described  by  Hill  [15]  in  connection  with  some  terts  run 
for  the  British  Broadcasting  Corporation. 


The  extension  of  this  mapping  to  higher  dimensions  can  be  accommo- 
dated simply  by  using  more  than  one  m-ary  pulse  to  subdivide  the  range  of 
m.  For  instance,  a bandspread -by-three  system  might  achieve  a stretch 
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factor  by  12  by  subdividing  the  range  of  m into  twelve  parts.  As  indicated 
in  Figure  9,  the  hg  pulse  could  be  3-ary  designating  which  third  of  the 
range  m occupies,  h^  could  be  4-ary  designating  which  fourth  of  the 
third  of  full  scale  m lies  in.  The  pulse  h0  would  then  carry  the 
remainder  information. 

If  the  range  of  values  of  hQ,  hv  and  ha , were  the  same  as  that  of  m 
ve  can  note  that  in  the  example  of  Figure  9: 


m = T7C4-3  -h, 


+ 4 


hi+  V 


12  h 
17  ^ 


4 , 

T7  hl  + 


where  the  transmitted  signal  is  the  n-tuple  [ha , hr  h0]  and  the  coefficients 
may  be  thought  of  as  the  basis  for  the  vector.  For  this  reason  the  tech- 
nique has  been  named  mixed  base  modulation.  The  selection  of  the  basis 
i s such  that  if: 


m 


B h + B , h , + 
n n n-1  n-1 


BP  h8 


+ Bj  hj  + B0  hQ 


(6) 


then: 

B^  = for  k 2 1 and  B = K 

where  K is  selected  such  that: 
n 

£ B.  = 1 
j = 0 J 

and  the  number  of  discrete  levels  of  h is  M , h , is  M ,,  etc.,  and  h„ 

n n n- 1 n-1  ° 

is  continuous.  It  is  interesting  to  note  if  through  are  equal  to  two 
the  system  is  equivalent  to  a binary  PCM  system  with  an  extra  bit  which 
describes  the  quantizing  error.  Such  a technique  was  suggested  as  early 
as  1948  by  Earp  [ 161  and  described  in  a paper  by  Nomura  and  Yasuda 
[171  in  1962. 


J 


14 


f 


s 

I 


4 


b 


in 

o 

n 


n 


RADIATION 


SECTION  V 


PERFORMANCE  OF  MBM  USING  A SIMPLE  RECEIVING  SYSTEM 


The  discussion  of  mixed  base  modulation  thus  far  has  considered 
only  the  modulator.  The  simplest  demodulator  is  that  shown  in  Figure 
10  where  each  noisy  discrete  pulse  is  quantized  to  select  the  nearest 
possible  transmitted  level  and  the  continuous  pulse  is  simple  left  alone. 
The  estimates  of  the  transmitted  values  are  then  appropriately  weighted 
to  provide  the  estimate  of  m.  This  is  equivalent  to  maximum  likelihood 
reception  except  for  an  end  effect  in  receiving  the  continuous  pulse. 


An  analysis  of  the  performance  in  terms  of  (SNR)  versus  (SNR),  is 
given  in  Appendix  2.  Performance  curves  are  given  in  Figures  11-14 
for  the  case  where  the  peak  power  in  all  pulses  is  made  equal  and  where 
the  number  of  levels  in  all  discrete  pulses  is  made  equal.  Bandspread 
factors  of  2 through  5 are  given.  The  probability  density  function  of  m(t) 
is  assumed  to  be  uniform.  The  rate  distortion  bounds  for  the  uniform 
source  are  obtained  from  equation  (3).  The  performance  of  a single- 
sideband system  is  plotted  also  for  reference.  Note  that  the  MBM  system 
trades  threshold  performance  for  high  signal-to-noise  ratio  gain  in  a 


manner  that  is  very  close  to  that  of  the  optimum  modulator-demodulator 
pair  which  allows  infinite  memory.  As  can  be  seen  from  Figure  11,  with 
a bandspreading  of  two,  40  dB  SNR  can  be  obtained  with  an  8-ary  and  a 


continuous  pulse  with  11  dB  SNR.  as  compared  with  40  dB  SNR.  for 


single- sideband.  The  ultimate  bound  is  between  2.2  and  Z3  dB. 


The  requirement  that  all  pulses  have  equal  energy  is  not  fundamental. 
The  division  of  energy  which  maximizes  SNR  for  a given  SNR.  is  determined 


o 1 

in  Appendix  4.  The  resulting  SNR  is  plotted  in  Figures  15  - 18.  As  can  be 


seen  in  most  cases  very  little  is  to  be  gained  by  unequal  energy  division. 
The  optimum  ratios  are  given  in  Appendix  4. 
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P4  OIATION 


PERFORMANCE  WITH  A MINIMUM- MEAN-SQUARE  ERROR  RECEIVER 


The  equation  for  the  minimum-mean- square  error  receiver  when 
bandspreading  is  present  (given  in  Appendix  1)  is  somewhat  cumbersome. 
For  the  case  of  MBM,  however,  it  is  shown  in  Appendix  3 that  the  mini- 
mum mean  square  estimate  of  the  transmitted  value  & given  the  received 
set  of  signals,  r,  - r , can  be  obtained  by  weighting  the  minimum  mean 
square  estimate  of  each  transmitted  pulse  given  the  corresponding  received 
signal.  This  argues  that  the  discrete  pulses  can  be  operated  upon  independ- 
ently and  if  the  discrete  pulses  in  the  transmission  have  the  same  number 
of  possible  levels,  the  same  operation  needs  to  be  performed  on  each 
pulse.  Of  more  significance  is  the  fact  that  if  the  receiver  minimizes 
mean- square  error  based  on  a uniform  density  function  for  the  input,  m, 
then  its  performance  will  be  better  than  that  of  the  simple  receiver  des- 
cribed in  the  preceding  section  for  all  density  functions  which  are  zero 
outside  the  range  of  the  uniform  variable.  This,  then  makes  the  minimum 
mean  square  error  receiver  of  practical  as  well  as  academic  interest. 

Unfortunately  the  calculated  performance  of  the  mmse  receiver  based 
upon  a uniform  density  function  (given  in  Appendix  4)  indicates  that  the 
performance  gain  is  small  (generally  less  than  a dB)  over  that  available 
from  the  simpler  receiver  described  in  the  preceding  section. 
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SECTION  VII 

COMPARISON  WITH  OTHER  TECHNIQUES 

Figures  19  through  23  are  performance  curves  for  SNR  versus 
SNR.  for  a variety  of  common  modulation  techniques.  In  all  cases  the 
same  idealizations  were  made  that  have  been  made  in  the  calculation  of 
the  MBM  performance  curves.  The  detailed  basis  for  the  curves  is 
given  in  Appendix  5.  A significant  point  to  note  is  that  in  all  of  these 
systems  the  high  SNR.  versus  improvement  threshold  tradeoff  is  tied 
directly  to  the  bandspread  factor,  g,  whereas  MBM  allows  the  trade  to 
be  made  for  a fixed  bandspread  factor. 

The  FM  performance  given  in  Figure  19  is  calculated  on  the  basis 
of  a peak-to-peak  frequency  deviation  equal  to  the  IF  bandwidth  and  a 
noise  model  based  upon  an  unmodulated  carrier.  The  values  of  g given 
are  the  ratios  of  the  IF  to  video  bandwidths.  Figure  20  shows  the  per- 
formance of  coherent  phase- shift  keying  and  quadraphase.  In  Figure  20 
the  values  of  SNR  for  large  values  of  SNR.  are  limited  by  quantization 
error.  In  Figure  21  three  different  coding-decoding  techniques  are  given 

which  provide  essentially  equivalent  bit-error  performance  versus  E /N 

b o 

( and  hence  equivalent  SNR  versus  SNR.)  but  requires  different  channel 
bandwidths.  The  values  for  g for  each  depend  upon  whether  PSK  or  QP 
are  used.  The  logic  in  listing  all  three  rather  than  the  most  conservative 
in  bandwidth  (the  63,  50  block  code)  is  that  the  23,  12  is  a more  commonly 
used  block  code  and  the  constraint  length  5 convolutional  code  allows  the 
simplest  decoding  equipment  of  the  three.  The  two  block  codes  are  BCH 
codes  and  a minimum  distance  decoder  is  assumed.  For  the  convolution- 
al code  of  Figures  21  and  22  the  decoding  strategy  is  assumed  to  be  that 
of  a Viterbi  decoder.  The  orthogonal  and  biorthogonal  systems  of  Figure 
23  are  based  upon  coherent  matched  filter  detection  and  hence  might 
apply  to  digitally  coded  systems,  discrete  PPM  or  banddividing  FM  such 
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as  described  by  Akima  [18]  (PFM  as  described  by  Rochelle  [19]  ) if  a 
coherent  reference  can  be  obtained. 

In  Figure  24  the  performance  of  the  various  systems  are  plotted  for 

a SNR  of  40  dB.  In  the  case  of  the  discrete  systems  only  the  minimum 

o 

value  of  g which  can  provide  40  dB  SNR^  is  shown.  Other  values  of  g 
are  either  larger  and  require  higher  SNR.  or  cannot  provide  40  dB  SNR^. 

Only  the  63,  50  block  code  from  Figure  18  is  shown.  As  can  be  seen  the 
bandspread  region  occuppied  by  MBM  is  not  available  from  the  common 
discrete  transmission  systems  and  the  continuous  systems  capable  of 
operating  in  that  region  are  not  very  efficient  in  their  use  of  these 
relatively  small  bandspread  factors.  As  the  bandspread  ratio  increases 
the  straight-forward  MBM  system  would  appear  to  be  • ’ikely  to  be 
competitive  in  performance  with  coded  quadraphase  but  might  still  be  a 
reasonable  choice  based  upon  equipment  consideration. 
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SECTION  VIII 

PRACTICAL  CONSIDERATIONS 

Several  practical  considerations  are  worthy  of  discussion.  The  MBM 

system  analysis  was  based  upon  use  of  cardinal  functions  for  pulses  with 

bandwidth  B 1 / (2 T ).  To  retain  this  bandwidth,  but  to  transform  this 

c o 

spectrum  to  a frequency  appropriate  for  propagation  it  is  necessary  for 
either  a single  sideband  technique  or  quadrature  modulation  (such  as 
that  used  for  the  color  subcarrier  in  color  TV)  be  used. 

If  one  takes  into  account  the  use  of  real  pulse  shapes  rather  than 

cardinal  functions,  the  bandwidth  required  in  the  receiver  preceeding  the 

matched  filter  may  be  much  wider  than  1/T  . This  implies  that  the 

spectral  occupancy  of  the  system  by  almost  any  definition  is  larger  than 

the  values  of  B given.  The  same  fact,  however,  applies  to  the  other 
c 

systems  discussed  in  Section  VII.  A thumb  rule  for  a PSK  system  which 
uses  square  modulating  pulses  is  that  a filter  of  4 times  the  bit  rate 
should  be  used  before  the  matched  filter  to  avoid  mismatch  and  crosstalk. 
The  spectral  occupany  of  FM  by  most  definitions  is  wider  than  the  IF 
bandwidth. 

It  is  reasonable  to  argue  that  the  crosstalk  problem  in  MBM  due  to 
bandlimitting  will  be  more  severe  than  that  of  a binary  transmission,  since 
MBM  may  use  pulses  with  more  levels.  For  that  reason  pulse  shaping  of 
MBM  pulses  will  almost  certainly  be  a requirement  and  must  be  consider- 
ed in  the  complexity  tradeoff  with  binary  (or  quadraphase)  systems  where 
pulse  shaping  is  rarely  employed.  However,  it  is  felt  that  reasonable 
pulse  shaping  techniques  can  be  employed  which  allow  use  of  IF  bandwidths 
and  provide  spectral  occupancies  which  are  closer  to  the  values  quoted 
than  is  the  case  for  quadraphase  and  bi-phase  systems  where  pulse  shaping 
is  not  employed.  Means  for  synchronization  of  the  receiver  demodulator 
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and  for  providing  timing  for  the  matched  filter  and  for  demultiplexing 
in  MBM  have  not  been  accounted  for.  Insertion  of  a reference  pulse 
periodically  could  serve  both  of  these  functions  plus  the  additional 
function  of  providing  a reference  for  calibration,  which  would  be  likely 
to  be  required  when  large  values  of  SNR^  are  desired.  In  general,  it 
is  felt  that  MBM  systems  can  be  constructed  such  that  the  total  imple- 
mentation losses  are  held  to  1 to  2 dB  which  is  a reasonable  estimate 
for  implementation  losses  in  the  other  systems  over  the  theoretical 
performance  curves  given  in  Section  VII. 
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DERIVATION  OF  THE  OPTIMUM  MEMORYLESS  OPERATORS 


Assume  the  two  operators  to  be  0T  (m)  and  0?  (r)  as  shown  in  Figure  2. 
The  optimum  choice  of  (r)  is: 

m = 0S  (r)  = E [m|  r] 

&=  f m f(m I r)  dm  = f m)  f<m)  dm 

> J f(r) 


& 


^ m f(r|m)  f(m)  dm 
^ f ( r | m)  f(m)  dm 


However,  f(r|  m)  = f [ r-Oj  (m)] 


-[r-01  (m)]2  /2a3 
e 1 n 

J2.v a 

n 


where  a2  is  the  variance  of  the  noise: 
n 


Thus,  the  optimum  receiver  for  any  0l  (m)  is  one  that  chooses  m to  be: 


A 

m 


° m -[r-a(m)]s  /2CT3  f 

\ m e n f(m)  dm 

F -rr-01  (m)]2  /2rr5 

\ e n 


= Xjr) 

D(r ) 


1-1 


For  such  a receiver  the  error  is  uncorrelated  with  the  estimate.  Thus: 


E f - m )3  = 0 

E [mrh]  = E [/h2] 


RADIATION 


and  if  E [ m2  ] = 1 

G = f (m-m)2  ] = E [m2]  - E[&2] 


Since  m is  dependent  upon  the  random  variable,  r,  the  expected  value  of 
m2  is : 

C m f(r|m)  f(m)  dm"!2  f(r)  dr 

E [ vr?  ] = C m2  (r)  f(r)  dr  = ^ 

> ’ f2  (r) 

1”  C m f(r|m)  f(m)  dm"!2 

E [Ai2]  = \ Jb dr 

^ ^ f(r  | m)  f(m)  dm 


Hence : 


e 


f(m)  dm  - 


r N2  (r) 
.)  D(r) 


dr 


1 -2 


By  use  of  standard  variational  techniques  under  the  constraint  that: 


P = ^ 0%  (m)  f(m)  dm 


1-3 


we  may  determine  that  the  optimum  choice  of  0X  (m)  satisfies: 


2 X (m) 


^0= 


r NUi 

\ D(r) 


[ r-0l  (m 


)]  T 2m 


Nilll 

D(r)  J 


-[  r-0.  (m)]2 /2cr2  dr  1-4 
e 1 n 
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whf  re 


r/  \ C r/  v -f  r-0.  (m)]2  /2a2  , 

I(r)  \ m f(m)  e 1 n dm 


and  \ is  a constant  which  is  chosen  to  satisfy  the  power  constraint  of 


1-3.  If  three  memoryless  operators  are  as  illustrated  in  Figure  3,  the 


optimum  operator,  03  (ra  , rp  ) is: 


^ m f ( rj  , r,  |m)  f(m)  dm 
f ( ri  • rs  ) 


r ([  r.  - 0.  (m)]3  + [ r_  - 0 (m)]3  )/2 as  ,,  . , N1,  . 

\ m e 1 1 L 3 3 n f(m)  dm  = N(rx  , ra  ) 


^ e([r,  - 0,  (m)P  + [r,  - 0,  (ml]/2o»n  f(m)  dm  ' D(Pl  , r>  , 


sing  the  same  procedure  as  before: 


c = Ff(m-m)3]  = E[m3]  - E[m3] 


G C rn3  f(m)  dm  - rj~"  \ \ ^ * drjdi 

\ ,/ZtTq  j .)  D(rj  , rs  ) 1 


This  procedure  easily  generalizes  to  N memoryless  operators  to  yield: 


m = 0 3 ( rt  , r8  , - - r ) 


N(rt  , ra  , - - rn) 


D( ri  . ra  . --rn) 


E 
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APPENDIX  2 


DERIVATION  OF  SNR  EQUATION  FOR  ONE  DISCRETE  PULSE 

o 


Let  us  analyze  the  output  signal-to-noise  ratio  of  an  MBM  system 

with  bandspread  by  two  and  mapping  as  shown  in  Figure  6.  Assume  that 

the  data  is  uniformly  distributed  between  _+  A and  that  this  range  is 

mapped  into  M lines.  The  two  dimensions  are  transmitted  as  two  pulses, 

the  first  with  uniform  distribution  over  a range  + a A and  width  T and 

— o o 

the  second  an  M-ary  pulse  with  M equally  likely  levels  ranging  from 
+ cLj  A and  width  T*  . 

The  receiver  is  assumed  to  make  a hard  decision  with  the  second 

received  pulse  as  to  the  line  corresponding  to  the  second  transmitted 

pulse  and  then  accept  the  first  pulse  directly.  If  the  received  value  of 

the  first  pulse  is  h and  the  discrete  decision  value  of  the  second  pulse 
o 

is  k,  then  the  estimated  transmitted  value,  m,  is 


& = k( 


2A>  r _o  ■ 

M J + a M 


Osk$  M-l 


If  the  actual  transmitted  value  was  m: 


. / 2A  \ 

= j ( ~rr  l + m - A 
\ M / o 


then  the  error,  e,  is: 


E • (m-Al  » (k-j»  ( ) ♦(^  - >"„) 


Howeve  r , 


h - a Mm  + n 
o o o o 
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where  n is  the  noise  associated  with  the  continuous  pulse, 
o 

Thus : 


_ . a M m + n 

„ . > / 2 A \ , f o o o \ 

€ = (k-jl  ( m“  ) * i ~ ■ r~M 

o 


o ; 


n 


€ = (k-j)  ( + ° 


M J a M 
o 


2-5 


The  two  error  terms  are  caused  by  the  additive  noise  on  the  two  channels 
and  are  independent.  Hence: 

o 

To  evaluate  the  variance  of  the  error  in  the  discrete  pulse  it  is  necessary 
to  evaluate  Ef  (k-j)3  1 


M M 

Ef(k-j)3]  =£  E (j-k)3  P(j|k)  P(k) 
j=l  k-1 


M M 

E[  (k-j)2  3 =—  E E (j-k)3  P(j|k) 
j=l  k= 1 


From  Figure  7 we  can  deduce: 


P(1 1 k) 


P(M|  k) 


f(|l-k|-i)2Aai  ^ 
Q<  °ni  (M- 1 ) I 

Q |(l  M-k|  -£)  ZAai 

' o " (M-l)  ^ 

n* 


j = i 


j = M 


2-7 
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(|  J ~ k I ~i )2Aa  ,(|j-k|+J)Aa 

p(jl  k)  = Q-  TV  \ - Q i 777-7-.  f 

l a (M-l)  J la  (M-l)J 

ni  m 


Z < j < M-l 


where 


00  - Z 


Q(x)  -\  ~ rr — dz  and  a is  the  standard  deviation  of  the 

1 JZ  tt  n 1 


noise  on  the  discrete  pulse. 


Since  E[(k-jP]  is  a function  of  the  absolute  value  of  j-k  only: 


M j-1 

E[  (k- j p ] = — T £ (j-kp  P(jjk) 
j=2  k=l 

We  have  dropped  the  terms  where  j = k. 


M-l  (M-k-%)2Aa 

E[k-jP]=  — E (M-kp  Q-j 1 j 

k=l  a (M-l) 

m 


M-l  j-1 


+ii  £2  = 

j=2  k = 1 L- 


0-k-|)2Aai 
~a  (M-l) 

ni 


(j -k+i )2Aa  ^ 

- Li  i 

■ m (M-l)  J * 


If  M-k  = m in  the  first  summation  and  j-k  = m in  the  second: 


Er(k-jP]  = E m2 

m=  1 


(m-i)2Aa 

Q | i\ 

l a (M-l)J 
m 


RADIATION 


M-  1 M-l-m 
+ — £ £ m2 

^ m=l  k=l 


[Q{ 


(m-^)2Aai  ^ 

(M-l)  J 
m 


- Q 


(m^)2AaM  , 

1 a (M-l ) j J 

XU 


M-l  (Zm-l)Aa 

EflMI3]  = m E Q{  ci  (M-l)} 

m=l 


m 


+ m z m ■ Qia  "(m^t;  J 

1 m = 1 - ni 


(2m-  l)Aa 


m=l 


M-2  (2m+l)Aa 

-ME,  (M-l  -m)  m2  I^Q  { “ (M-1)J  J 


m=l 


m 


Letting  m = n-1  the  last  summation  and: 


(2m-l)Aai 

Q{~0  (M-l)J  ~ °m 

XU 


M-l  M-2 

E[(k-jf  ] = — {E  m2Qm  + £ (M-l-m)m2  Q] 

^m=l  m=l 

M-l 

-£  (M-n)(n-l  f Q 
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= h { Q>  {U-lfQu-i  * - <m-2|2qm-i 

M-2 

+ 2 [m2  + (M- 1 -m)3  - (M-m)(m-l  J3”] Q /• 

m = 2 -«  n/ 

2 M'2 
= M {(M-1)Q1  + (2M_3)Qm-1  +S  (2m-l)(M-m)Q  1 

m=2 


Ef(k-jf] 


2_ 

M 


M-l 

S 


m = l 


(2m- 1 )Aa 

(2m-l)(N-m)Q{— 

ni 


Thus  the  mean  square  error  is: 


E[e2 ] = 


no 

a2  M2 
o 


8 A2 
M3 


M-l 


(2m- 1 )Aa 


(2m-l  )(M-m)  Q 


m=  1 


l a (M-l ) 

ni 


4 


The  signal  component  of  (ri  is  equal  to  m.  Hence  the  mean  square 
output  signal  to  noise  ratio  is: 


n 


A^ 

3E[e] 


1 

M-l  (2m- 1 )Aa 

3a  3 /a2  M2  A2  t24/MBi;  (2m- 1 )(M-m)  q{ 

no  o ' ' lo  (M-l) 


The  power  in  the  first  pulse  is: 


P 

o 


a2  A2 
o 

3 


and  the  energy  is: 

a A3  T 

E - ° ° 

o 


2-10 


2-11 


3 
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The  variance  of  the  noise  adding  to  the  first  pulse  is: 


N 

j 2 = — — 

no  2T 


The  variance  of  the  second  pulse  is: 


M-l  2Aa 


4 A2  a2 


P = 1 — 

1 M(M-l)  (_ 


2 A2  a2  (2M-1! 
Pl  “ 3(M- 1 ) 


tjr  r <M=UM^M-JLT|  . (Aai  f 


p _ _l f M + In. 

1 3 V M-l  ) 


The  ene  rgy  is  : 

E = P T = JL—L  ( M.+JA 
1 11  3 V M - 1 ) 


The  variance  of  the  noise  in  the  second  pulse  is: 


m 2T. 


The  energy  in  the  total  signal  is. 


a?  A2  T w . , a2  A3  T 

E - E + E — - ( MJLL\  + _2 

Es  " Eo  El  ‘ 3 Im-1  3 


E = — Ta3  T.  ( + a3  T *1 

s 3 1 \ M - l ) o oj 
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If  we  define  : 


E 
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P‘  ~f>  • tf  >*a««  %!fP 


a2  T 
o o 


° E P ,M+1 , p 

8 a?  T,  — ) + a2  T 
1 M - 1 o o 


E.  4T,<|£rf> 

Es  af  Tx  $±j-)  + a2  T 
1 M- 1 o o 


then: 


l/M2  (—“)  + 

o 


24/M3  1 (2m - 1 )(M-m)Q|  (2m-  1 )'i~  jf-  (^~ 


if  2 E /N  = Y 
s o 


■’“n  1 + 5±5ax 

M 


M2  RoY 


(2m-l  )(M-m)Q|(2m-l  )(3R1  Y/M2  -1)*} 


The  case  where  multiple  discrete  pulses  are  utilized  may  be 

analyzed  in  a similar  manner.  Let  us  assume  that  J discrete  pulses  are 

transmitted  with  M.  levels  for  the  such  pulse.  The  peak  values  of  the 
.th 

j pulse  are  _+  Aa^.  As  before  the  peak  values  of  the  continuous  pulse  is 
Aa^.  The  first  discrete  pulse  is  used  to  divide  the  total  range  into 
•egments ; the  second  divides  each  of  the  segments  into  smaller 
segments,  etc.  If,  as  before,  decisions  are  made  in  the  usual  manner 
on  the  discrete  pulses,  the  error  at  the  receiver  is: 


t 


n 
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e = 


a u M. 
o J = 1 j 


, 2A_  . . 2A 2A 

M,  ,k>  1 M,  M,  - 


r M. 
J=1  J 


'j!5 


e = 


+ £ 


Z A 


a tr  M.  n=l  t/  M. 

° J=1  J J=1  J 


(k  -j  ) 
n n 


Z-18 


where  as  before  the  level  of  the  pulse  was  transmitted  and  the  k^  level 
received.  The  mean  square  error  is: 


E[e3]  = 


no 


asOJ  ,M.]: 

o Lj  = l jJ 


J E[(k  -j  ry 

+ 4AS  2 

n=l  (v  . M.)2 
J = 1 J 


2-19 


Using  equation  2-10  to  express  the  ET(k  -j  f]  we  have: 

n n 


a J Mn-1 

no  n 

31  = r + 8A3  S 1/M  (t  ,M.f  2 (2m-l  )(M 

2 r nip  i n 1 1 J i ^ 

a2[ir  ,M.r  n=l  J J m-1 

o J = 1 J 


E[€2]  = 


-m) 


(2m- 1 )Aan 
(M-1 ) 


Q{^ 

La 

nj  n 
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The  signal  to  noise  output  is  then: 


o 

T2 

n 


J J n Mn-1 

3cno3 /Aaa03  ^jT=1Mj-'3  +24^  (l/Mn(tr=1M.)a)E  (2m-l)(Mn-m) 

n = l m = l 


where: 


Q = Q ji_2 
mn  L a . 

nj  n 


(2m -1  )(Aan  l 

J 
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APPENDIX  3 

PERFORMANCE  OF  THE  MINIMUM  MEAN  SQUARE  ERROR 
RECEIVER  FOR  MBM 

The  equation  for  the  MMSE  receiver  for  the  general  N + 1 
dimensional  memoryless  mapping  given  in  Appendix  1 is  rather  cumber- 
some. The  MBM  form  of  mapping,  however,  reduces  the  problem 
considerably.  This  can  be  argued  in  the  simplest  fashion  by  noting  that 
the  input  is  in  the  form  of: 

N 

m = H B.  h 3-1 

i = 0 1 1 

The  received  signal  set,  r , r r , are  related  to  the  trans- 

o n 

mitted  signals,  h , h, h.,,  in  that: 

o 1 N 

r.  = h.  + n. 

J J J 

The  random  variables  h.  and  n.  are  jointly  independent.  Thus: 

f(h  , \ , h , r , rj r ) = f(h  ,r  ) ffl^  , r*  ) f(h  ,r  ) 

o 1 No1  n o o ^ 1 NN 

and  hence : 

f(ho,  h,  , hN|  rQ,  r*  , rN)  = f(hj  tq)  f(h,  | r,  ) f(hN|  rN>  3-2 

The  optimum  receiver  is  one  that  outputs: 

& = E[m|  r , ^ , r ) = \ m f(m|  r , rt rN)dm 


H 


PfA  D/A  T/ON 


™ Bihi)f(ho*N— -hNlro*  ri—  -rN)dhodhi—  -^r 

N 

& = (S0  Bihi)  f(hol  ro)  f(h^  I r* ,— -f(hNl  rN)  dho  dh> dh; 


which  yields: 


IN  IN  AN 

m = £ B.  Ch  f(h.|r.)dn.  =E  B.  E{h.  I r.  ) = £ B.ti.  3- 

. _ 1 ) i i'  i i . A l l i i 

i=0  J i=0  i=0 

The  equation,  3-3,  states  that  the  minimum  mean  square  estimate  of 
m is  a weighted  sum  of  the  minimum  mean-square  estimate  of  each 
transmitted  pulse.  Further,  the  error  in  the  estimate  is: 


m - fn  = £ B.  (h.  -ti. ) 

i=0  1 1 1 


and  the  mean  squared  error  is: 


N N 


e E(m2  ) = E[  (m-m  f ] =£  £ B.  B.  E-f  (h.  -ti. ) ( h.-Ji.A 

i=0  j=0  1 J 1 1 1 J J f 


Since : 


E[h.h  1 = 0 and  E [h.ti.]  =0  i t j 

i j 1 J 

N N 

e EJm2  ) = E[(m-mf]  = £ B2  E^(hj-?i  f j = £ B2  e2 


Thus,  the  mean  square  error  associated  with  the  MMSE  receiver  is  the 
weighted  sum  of  the  mean  square  error  associated  with  the  optimum 
estimate  of  the  transmitted  values. 


A . Hit) 


E 
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tf  AM.  * :AMH  S A >H 


Ln  the*  case  of  the  continuous  pulse: 


f(V  = zTa 

o 


-a  A < h.  < + a A 
o 1 o 


From  Appendix  1 


+a  A 
^ ° h 

-a  A 


-(r  -h  )2  IZn  2 ,, 

e o o n an 

o o 


h = C h f(h  I r ) dh  = — -2- 

o ' o o'  o o +a  A , 

r o -(r  -h  )2 /2a  2 
\ e o o n 


-a  A 
o 


Performing  the  integration  yields: 


a ^ -(a  A-r  )2 /2 a 3 -(a  A+r  f /2a  21 

ti  = r - — -co  ° n-e  o o nj 

o o ./2rr  /a  A-r  \ /-a  A-r  \ 

Qi  o o ) - Q ^ o o ; 

a ' a 


*nce:  r = h + n 
o o o 


li  = h + n - -7== 
000  J Ztt 


0 j”  -(a  A-h  -n  )a /2 a 2 -(a  A+h  +n  )2 /2a  2 

nleo  00  n-eo  00  n 


_/a  A-h  -n  \ ^/-a  A-h  -n  \ 

Qi  o 00  j - Q ( o 00) 


n = h + n - F(h  , r ) 
000  00 


The  mean  square  error,  e 2,  is: 

o 


e 3 = E[(h  -h  )21  = E[n  2 1 - 2E[n  F(h  ,n  )]  + E[  F®  (h  ,n  )] 
o 00  oJ  ‘•o  00  L 00 


e 3 = o 2 - r^-r  K f 2n  F(h  , n ) - F®  (h  , n )]  °n — - 
o n 2 A \ o'oo  o oJ  ./Tna 


n 2 / 2 

0 ' dh  dn  3-6 
= o o 
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where : 


F(h  ,n  ) = 75= 
o o v2u 


r -(a  A-h  -n  )3 /2o  2 -(a  A-h  +n  )3 /2 a 2_j 

* - a o o nj 


n -e 


Q 


/ a A-h  -n  \ ^ / -a  A-h  -n  \ 

(-2 2)  - Q( — ° 2_e) 


n 


Although  in  equation  3-6  it  is  desired  to  average  over  the  value  of  h 

o 

as  shown,  it  is  interesting  to  observe  the  value  of  e 3 for  a fixed  value 

o 

of  h . If  the  ratio  of  peak  signal  to  rms  noise,  a A/a  , is  reasonably 
o on 

large  (2  5)  then  the  calculatable  ratio  of  e 3 to  a 2 depends  only  upon 

on 

the  number  of  noise  variance  intervals  that  h differs  from  plus  or  minus 

o 

a A.  The  reduction  in  e 3 is  plotted  as  a function  of  a 3 in  Figure  3-1. 
o o n 

The  MMSE  receiver  based  on  a uniform  density  does  not  increase  the 

error  for  any  transmitted  value  of  h . This  fact  can  be  verified  from 

o 

equation  3-4.  Therefore,  such  a MMSE  receiver  will  provide  results 

with  the  continuous  pulse  which  are  at  least  as  good  as  those  provided 

by  the  simpler  receiver  for  any  density  function  which  is  zero  outside 

the  range  a A and  -a  A.  The  average  reduction  in  error  for  values  of 
o o 

h such  that  la  A-h  I £ 5a  is  0.181  a 2 • Since  the  reduction  when 
o ' o o'  n n 

hQ  is  outside  this  region  is  negligible  we  may  deduce  that: 


5a  _ 


= aj 3 [1  - .181  (—5)] 


n 


a 2 A2 

If  Y = R Y = . v a 
o o 12a 


a A' 
o 


a A 2 8a 
o n 


then:  e 3 =a  2 [1  - . 905/(12Y  )*]  = a 2 [1  - .262/Y^l  Y 2 5.33(7.3dB) 
onL  o n 0 0 

Figure  3-2  is  a plot  of  e 3 versus  Y for  a MMSE  receiver  and  the  simple 

00 

receiver  described  in  Section  V. 
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In  the  case  of  the  discrete  pulse,  h,  with  M levels  ranging  from  plus  to 
minus  A,  the  MMSE  receiver  is: 


M a A 

h(r)  = E (h | r)  = £ h P(h.|  r)  where  h.  = -77 — ( j - 7 ^ 

j = 1 J J M \ 2 / 


S.  h,  P(h.|  r) 
i = l 1 r 


£,  f(r|h.)  P(h.) 
J = 1 J J 


Since : 


-(r-h.)3  /2CT  3 

f (r  j h. ) = -==—  J 

J V 2ir  a 


P(h)  = r/ 

j M 


then: 


^ , -(r-h.)3  /2a  3 

a h e J n 

h(r)  = L- }- 

^ -(r-h.)3 /2a  3 

ple  J 


If  the  actual  transmitted  value  was  h then  r = h + n and: 


h(r)  = 


™ -[(h-h.)  +n]3/2a  3 

u , h . e k 1 n 

1 = 1 1 J 


e-[(hk-h  ) + n 

j = l 


]3/2a 
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If  the  k term  in  each  summation  is  taken  out  of  equation  3-6 


-n2  /2o  3 T -[(h  -h.)  + n]2  /2a  2 

h,  e n+Zj  h.  e ki  n 

k 1 = 1 iV  l J 


•n2 /2a  3 


3 -tOv-kO  +n]3/2o  3 

” +j5wk  e V J 


This  is  in  the  form: 


a a+b  a , b-ad/c 

?>(p)  = ^ = r + 


Thus : 


E,  (h.-hlr)e-[(hk-hj)  +n]3/2an3 

V)  = h + 

E e^-V+n^V 

j = I 


ti(r)  = h + F(k,  n) 
k 


The  second  term  is  obviously  the  error.  The  mean  square  value  of  this 
error  if  all  values  of  k are  equally  likely  is: 


1 M 

= M El  \ f*<k.n)d„ 


ra  A (j-k)  + n-|2  /2a  3 

I ~TT  n 


where  F(k,n)  = 


(~5r)  <j-k> 3 M 

M ra  A (j-k)  + n"13  /2a  3 

j n 


Equation  3-7  was  numerically  integrated  to  produce  curves  of  cp 


versus  Yi  for  values  of  M from  1 to  10  where: 


4 


r 


a 


Y -iL  JL 

Yl  " 3o  2 

n 

These  are  plotted  in  Figure  3-3  along  with  the  equivalent  error  produced 
by  the  simple  receiver.  As  was  the  case  with  the  continued  pulse  error 
evaluation  the  error  for  the  MMSE  receiver  given  a particular  trans- 
mitted value,  h^,  is  equal  to  or  less  than  the  error  with  the  simple 
receiver  for  all  values  of  k. 

Equation  3-9  can  be  used  in  combination  with  the  data  in  Figures 

3-2  and  3-3  to  calculate  SNR.  versus  SNR  for  MBM  with  a MMSE 

i o 

receiver  and  a specified  energy  ratio  for  the  pulses. 


■I 


n 

4 i 
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OPTIMUM  ENERGY  RATIOS 


The  mean  square  error  for  MBM  transmission  can  be  expressed 
in  the  form: 


E B2  e.3  (Y.) 
1 1 


i = 0 


4-1 


where:  J + 1 is  the  bandspread  factor,  g,  or  the  number  of  output 
pulses  per  input  pulse. 


B.  is  the  weighting  or  basis  for  the  ith  pulse 


2E. 


Y.  = 


i N 


= average  pulse  energy  to  noise  spectral  height  ratio 


for  the  ith  pulse. 


Equation  4-1  applies  to  both  the  case  of  the  simple  receiver  and  the 
MMSE  receiver.  If  the  total  energy  to  noise  ratio  is: 


2E 


Y = 


N 


= £n  Y- 
i = 0 i 


4-2 


then  the  optimum  problem  is  one  of  minimizing  e3  in  equation  4-1  under 
the  constraint  of  equation  4-2.  Fortunately,  this  problem  can  be  handled 
as  a sequence  of  two-dimensional  optimization  problems  rather  than  a 
J + 1 dimensional  problem.  If  one  starts  with  the  first  two  terms: 


e 3 = B 3 e 3 (Y  ) + B^  e 3 (Y,  ) 
mi  o o o 11  1 


PRBCroiNS  PAOS  BLANK-NCT  PTUED 
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and  determines  the  minimum  achievable  e 2 for  a given  Y,  ' = Y + Y,  . 

mi  1 o 1 

then  a new  curve  exists  such  that: 


e 2 = e,  ' (Y,  ' ) 


If  one  then  adds  the  next  term: 


= "mi"  + V ^ (Y=  ) = el  ' (Y,  ' ) + B,2  es2  (Y,  ) 
and  minimizes  for  all  Ys  ' - Ya  + Yx  ' this  produces  a new  curve: 


e 2 = ea'  (Ys  ' ) 

mp  •*  ** 

which  corresponds  to  the  optimum  split  between  the  first  three  pulses. 
This  procedure  can  be  continued  until  all  J + 1 pulses  are  included. 

Plots  of  ej2  (YJ  are  given  for  both  the  simple  and  MMSE  receiver 
in  Figures  3-2  and  3-3  in  the  previous  Appendix.  Since  the  analytical 
forms  for  these  plots  were  cumbersome  it  was  found  easiest  to  produce 
the  combined  optimums  by  hand  and  replot  the  new  curve.  The  resulting 
optimum  energy  ratios  are  given  in  Figures  4-1  through  4-4.  The 
resulting  SNR^  versus  SNR^  for  optimum  energy  split  are  plotted  in 
Section  IV. 
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SYSTEM  PERFORMANCE 


SINGLE  SIDEBAND  (SS) 

A single  sideband  system  transmits  data  of  bandwidth,  B . The 

v 

sampling  rate  is  2B  and  sample  time,  T is: 

V s 


TS  " 2B 


If  the  pulse  has  average  power,  P , the  energy,  Eg  is: 


E = P T = — — ■ 
s s s 2B 


If  the  pulse  is  transmitted  in  the  presence  of  white  noise  with  two  sided 
height  of  N^/2,  and  a matched  filter  with  impulse  response: 


H(t)  - 


0 < t £ T 


otherwise 


is  used  then  the  output  due  to  signal  is  equal  to  the  input  and  has 

variance,  a 3 > 
o 


a 3 = P 
o s 


The  variance  of  the  noise,  o 2 is: 

n 

NT,  N 

a = C s (-L_)8  dt=-2 

n 2 \ ' T ' 2T 

s s 

o 


Thus,  the  output  signal  to  noise  ratio  is: 


« to  - *.  \ % 
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SNR  = CT  2 /ct  3 
o o n 


2P  T 
s s 


N 

o 


2E 


s 


N 

o 


P 

s 

N B 

O V 


SNR. 

1 


5-5 


Since  the  output  sampling  rate  equals  the  input  sampling  rate  the  band- 
spreading  is: 


P = 1 


OTHER  FORMS  OF  AMPLITUDE  MODULATION 

Suppressed  carrier  A.  M.  differs  from  single  - sideband  only  in  the 
method  of  frequency  translation.  The  effect  is  to  leave  the  relationship 
between  SNR^  'aid  SNR.  the  same  but  double  the  bandspread  ratio.  That 
is  for  SCAM: 


SCAM  I SNR  = SNR. 

' o 1 

P = 2 

Ordinary  non-suppressed  carrier  A.M.  is  the  same  as  SCAM  except 
that  additional  power  is  added  to  the  signal  in  the  form  of  carrier.  This 
signal  increases  the  required  input  power  without  increasing  the  SNR^. 

For  the  case  of  a uniformly  distributed  modulating  signal  and  100% 
modulation  index: 

AM  I SNR  = .25  SNR. 

1 o 1 

P = 2 

FREQUENCY  MODULATION  (FM) 

The  general  equation  used  to  calculate  SNR^  versus  SNR.  for  frequency 
modulation  is: 
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1 2 o fd2  B ,/B  (12B  2 + B 2 ) 

if  V S V 

SNRo  ' 1 + o F(p)  • 3B  (B  -B  )/ (1 2B  2 + B 3 

if  if  v s v 


where : 


p = — — SNR.  = IF  signal  to  noise  ratio 

* c i 
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fd  = rms  frequency  deviation 

B = bandwidth  of  video  filter 
v 

Bg  = center  frequency  of  video  filter 

B.  = bandwidth  of  I.  F.  filter 
if 


F(p  ) = {(l+«)  Q (v/2d  ( 1 +z ) ) Q(-V/Tp(l  +z)) 


+ [z(l+z)]£  e D[Q(-yZpz)  - Qf+^Toz)]  [Q(-y?5(l +z)-Q(y2o(l +z)] 


+ ze  QC/Zpz)  Q(-x/?pz) 


whe  re : 


QM--{  b 


CD  -t2  /2 


f = explained  in  the  next  paragraph, 
dc 

Equation  5-6  is  taken  from  reference  [20]  where  the  term  represent- 
ing the  spectral  height  of  the  FM  output  noise  at  the  frequency  origin  has 
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been  replaced  by  F(p).  The  expression  for  F(p)  is  taken  from  reference 
[21]  and  is  the  spectral  height  of  the  noise  at  the  origin  based  on  the  noise 
model  of  that  report.  This  spectral  height  is  derived  assuming  an  unmodu- 
lated but  not  necessarily  tuned  carrier.  The  term,  f , represents  the 

dc 

detuning  of  the  carrier  from  the  center  of  the  IF. 

Equation  7-1  applies  to  cases  where  a bandpass  video  filter  is  used 
as  well  as  the  case  where  a lowpass  filter  is  used.  Figure  5-1  is  a plot 
of  normalized  SNRq  versus  SNR^  with  thresholding  noise  based  upon  no 
frequency  offset.  These  curves  for  the  case  of  the  lowpass  video  filter 
correspond  very  closely  to  curves  based  on  Stumpers  [22]  equations 
(see  Akima  [18]).  Figure  5-2  is  similar  except  that  the  threshold  noise 
is  based  upon  a frequency  offset  of  half  the  IF  bandwidth. 


Figure  19  in  Section  VII  is  taken  from  Figure  5-1  where  f^  = B 

and  = B^/ 2.  For  uniformly  distributed,  m,  this  corresponds  to  setting 

the  peak-to-peak  deviation  equal  to  the  IF  bandwidth.  For  FM  we  call 

B /B  =g. 
if  v 


DIGITAL  TECHNIQUES 
Phase  Shift  Keying  (PSK) 

The  bit  error  probability,  p , for  PSK  is  well  known  to  be: 

e 


2Et 

Pe  = Q C FT 


where: 

E,  = energy  per  bit 
b 

N^/2  = noise  spectral  density 


RADIATION 


o 


If  the  message  value,  m,  is  quantized  to  n bits  then: 


2E  2nE, 

s b 

SNR.  = — = — 

1 N N 

o o 


SNR. 


p =Q  C— “1 

e n 

Also  if  the  distribution  of  m is  uniform,  the  quantization  error  normalized 
to  the  rms  value  of  the  data  will  be: 


E 3 =1/2 

q 


n 


The  normalized  link  error  is: 


V = 1 - F"  ] 
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SNR 


o E 3 + E 


t l/2n  + 4ri-l/23ni  Qr^^i] 

j l n J 
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The  pulse  shapes  can  in  theory  be  made  so  that  the  channel  bandwidth 
is  equal  to  the  bit  rate  which  is  n times  the  sample  rates  and  2n  times 


B . 

v 


p = 2n 


Quadraphase  (QP) 

The  only  difference  between  quadraphase  and  PSK  is  in  the  channel 
bandwidth  which  is  one-half  that  of  PSK.  Thus: 

8 = n 

and  SNRq  versus  SNR^  is  given  by  equation  5-5. 


»s  c • 
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Coded  PSK  Systems 

The  calculation  of  bit  error  rates  as  a function  of  E,  /N  for  several 

b o 

error  correction  codes  is  shown  in  Figure  5-3.  As  before  the  relation- 
ship between  E./N  and  SNR.  is: 

bo  1 


2nE 

SNR.  = 5-9 

i N 

o 

where  2n  is  the  number  of  quantization  increments  of  the  code.  A fact 
not  widely  known  is  that  the  relationship  between  bit  error  probability 
and  rms  link  error  for  most  coded-systems  is  the  same  as  that  given  in 
equation  7-4.  The  necessary  and  sufficient  requirements  for  this  to  be 
true  are: 

(a)  The  probability  that  a particular  information  bit  is  in  error 
is  the  same  as  that  of  any  other  bit  being  in  error. 

(b)  The  error  in  a particular  bit  is  equally  likely  to  be  in  either 
polarity. 

The  argument  to  justify  the  statement  is  to  observe  that  the  mean 
square  error  when  several  bits  in  a single  word  are  in  error  is  the  same 
as  that  when  the  bits  occur  in  separate  words.  This  is  a consequence  of 
(b).  If  the  absolute  error  in  a particular  word  is  e j and  the  error  due 
to  an  additional  bit  error  in  the  same  word  is  which  is  equally  probable 
to  add  or  subtract  from  et  then  the  mean  squared  error  is: 

E(ea)  = . 5(e!  + % f + . 5^  - f = ex3  + e3a  5-10 

This  is  the  same  squared  error  that  would  occur  if  the  additional  error 
had  been  in  a separate  word.  If  all  bits  in  the  n bit  word  are  in  error 
with  probability  p and  the  mean  squared  error  is  the  same  as  that  which 
would  occur  if  they  were  in  separate  words  then  the  total  mean  squared 
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error  is : 


mse  = p 


n 

E 

j=i 


I* 


[1  - l/22n] 
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This  error  is  normalized  to  full  scale.  Normalizing  to  the  rms  value 
of  a uniform  density  required  multiplication  by  12  which  produces 
equation  5-8. 

Although  both  block  and  convolution  codes  tend  to  cause  clustered 
errors,  the  usual  decoding  techniques  still  meet  requir  mert  (a)  even 
in  the  case  where  blocks  are  an  integer  number  of  word 

Thus,  the  equation  for  SNRq  for  coded  PSK  is: 


SNRo  1 /23n  + 4 p (SNR.) 

e l 

A plot  of  p (E/N  ) is  given  in  Figure  5-3  for  several  coding -decoding 
e b o 

techniques. 

The  bandspreading  is  the  bandspreading  caused  by  the  coding,  a, 
times  the  bandspreading  of  the  PSK  system.  Thus: 


g = 2an  for  coded  PSK 

8 = an  for  coded  QP 
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Orthogonal  and  Biorthogonal  Systems 

The  symbol  error  probability ,pg,  as  a function  of  E^/N^  for 
orthogonal  and  biorthogonal  systems  is  well  known  and  has  been  tabu- 
lated by  Viterbi  [23].  Totty  [24]  has  shown  that  the  mean  square  error 
normalized  to  full  scale  of  orthogonal  signals  is: 


5-14 


E2 
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mse  = P /6 
s 

and  that  for  biorthogonal  signals  the  mean  square  error  is  very  nearly 
the  same.  Actually  the  argument  associating  bit  error  probability  with 
mean  square  error  given  in  the  preceding  section  can  also  be  applied  to 
orthogonal  and  biorthogonal  systems.  This  would  allow  exact  calcu- 
lation bor  biorthogonal  systems  since  the  bit  error  probability  as  well 
as  the  symbol  error  probability  was  tabulated  by  Viterbi.  Unfortunately 
the  tabulated  bit  error  probabilities  for  the  biorthog‘~al  case  are  in 
error  so  that  we  must  settle  for  use  of  equation  5-14  and  Totty's 
bounding  arguments.  Normalizing  equation  5-14  to  the  mean  square 
output  for  uniform  data  we  have: 


Orthogonal  and 
Biorthogonal 


o~  l/23n  +2  P (SNR.) 

s 1 


where  as  before: 


2 = number  of  quantizing  increments 

SNR.  = 2nE  /N 
i bo 

The  bandspread  factor  is: 

Orthogonal  g = 2° 

Biorthogonal  g = 2° 


■i 


FIGURE  3-2  ERROR  IN  CONTINUOUS  PULSE 
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